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Abstract 

In this study, a virus strain designated as HYl 2 was isolated from cattle with a disease of high morbidity and mortality in Jilin 
province. Biological and physiochemical properties showed that HYl 2 isolates is cytopathic with an extremely high 
infectivity. HYl 2 is resistant to treatment of organic solvent and acid, and unstable at 60°C for 1 h. Electron microscopy 
observation revealed the virus is an approximately 22-28 nm in diameter. The complete genome sequence of HYl 2 consists 
of 7416 nucleotides, with a typical picornavirus genome organization including a 5'-untranslated region (UTR), a large single 
ORF encoding a polyprotein of 2176 amino acids, and a 3'-UTR. Phylogenetic analysis clustered HY12 isolates to a new 
serotype/genotype within the clade of enterovirus E (formerly BEV-A). Alignment analysis revealed a unique insertion of 2 
amino acid residues (NF) at the C-terminal of VPl protein between aa 825 and 826, and several rare mutations in VPl and 
VP4 of HYl 2 isolates in relation to known bovine enterovirus (BEV) strains. This is the first report of an enterovirus E in China, 
which is potentially associated with an outbreak in cattle with severe respiratory and enteric diseases. 
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Introduction 

The genus enterovirus within the family picornaviridae consists of 9 
species of enterovirus (A,B, C, D, E, F, G, H, J) and 3 species of 
rhinoviruses (A, B, C) based on tlie latest virus taxonomy[l] . These 
viruses have many features in common and are the leading 
etiological agents related to respiratory and digestive diseases in 
human and animals. Like other species within the genus, bovine 
enteroviruses (BEVs) are small, non-enveloped viruses with an 
icosahedral virion and a positive-stranded RNA genome. BEVs 
have been isolated from cattle with a clinical signs varying from 
respiratory diseases to enteric, reproductive disease and infertility 
[2,3,4,5,6,7,8,9,10,11,12,13,14], and even from faeces of the 
presumably healthy calves [15]. The pathogenicity and virulence 
of bovine enterovirus is still largely unknown. Failure to 
experimentally reproduce bovine enterovirus infection in calves 
showing obvious clinical signs led to the conclusion that BEV is not 
significant agent in the cattle industry. This argument was seemly 
further supported by the findings that bovine enterovirus were 
easily detected in the contaminated waters adjacent to cattle herds; 
and the discovery that BEV-like sequences are present in shellfish, 
bottlenose dolphins, and in deer feces from the same geographical 
area [16,17,18,19,20]. However, with more and more BEV 
isolates were identified from the fatal enteric and respiratory 
disease, the pathogenicity and virulence of BEV recurred where its 
relatedness to the illness be intensively explored [12,21,22,23,24]. 
Blas-Machdo has showed that calves experimentally infected with 



BEV- 1 manifested symptoms of respiratory illness similar to those 
naturally BEV-infected calves; and the BEV-1 was detected to 
localize mainly in the digestive tract, indicating the potential 
pathogenicity of BEV [22]. Recendy, several strains of BEV were 
isolated from cattle with a severe diarrhea in China, and genomic 
sequence analysis indicated that those strains belong to enterovirus 
F (previously named BEV-B)[12]. Those findings suggest the 
pathogenicity of BEV as potential causative agents for the 
respiratory and enteric diseases in cattle. 

Classification of bovine enteroviruses undergoes a series of 
modifications. In the early attempt, BEV were classified into seven 
serotypes, and later revised to two serotypes [25,26]. Because of 
the cross-reaction among BEV type-specific sera, it is difficult to 
type BEV using serological means. With the accumulation of BEV 
sequence data, classification of BEVs based on virus genetic 
variability and molecular difference becomes feasible. Based on 
the generally accepted definition for picornavirus species and 
serotype, a molecular-based BEV classification by comparing the 
sequences from 5' -UTR, the capsid protein region were estab- 
lished and used to classify bovine enteroviruses to BEV-A 
(currendy named enterovirus E) and BEV-B (currently named 
enterovirus F), where different serotypes/genotypes were further 
divided for either enterovirus E or enterovirus F [27,28,29]. In this 
study, we reported the identification of a novel enterovirus E 
isolates HYl 2 from a cattle herd with an outbreak of a severe 
respiratory disease and enteritis in JUin Province. 
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Table 1. Primers used for differentiation of the potential agents. 





virus 


Primer sequences (5'-3') 


Position 


BPV 


S GCGCCGCATAAATGTGTCTTGGTG 


1326-1349 




AS TGTGGGGCTTTCCGCTTTATCTCA 


1738-171S 


FMDV 


S CGGTGGAAAACTACGGGGGAGAG 


40-63 




AS TGCGCCGTAGTTGAAGGAGGTTG 


496-473 


BEV 


S AGGATGATGATTGGCAGATTTTGT 


372-395 




AS CATGTGGAAGTGTCTTTTGAGGAA 


708-685 



BPV: bovine parvovirus; FMDV: foot and mouth disease virus; BEV: bovine enterovirus; S: sense; AS: antisense. 
doi:1 0.1 371 /journal.pone.0097730.t001 



Materials and Methods 

The ethics committee in JUin University has approved this 
study. 

Cell culture and virus isolation 

Faecal samples were collected from diarrheic cattle following 
standard procedures approved by the ethics committee of College 
of Veterinary Medicine at Jilin University, and processed to 
inoculate the MDBK cells. Briefly, samples were diluted in a 
dilution of 1:10 (W/V) with 10 mM phosphate buffered saline 
(PBS) (pH 7.2). After centrifugation at 8000 r. p. m. for 30 min, 
the supernatant was filtered with 0.45 nm filter and the flow- 
through was used to inoculate MDBK cells. The inoculum was 
discarded after incubation with MDBK cells for 2 h, and the cells 
were washed with Hank's solution before the addition of 
Dulbecco's modified eagle's medium (DMEM) (Invitrogen,) 
supplemented with 2% fetal bovine serum (HyClone), 2 |J.g/ml 
gentamycin and 2 mM L-glutamine (Invitrogen). The cells were 
observed every 4~6 h and cytopathic effects (CPE) were captured 
using Digital Camera. 

Electron microscopy observation 

Sample was processed for EM observation by centrifuging at 
8000 r. p. m for 30 min at 4°C after virus cultures were frozen and 
thawed for 3 times. The supernatants were incubated with 1 % 



phosphotungistic acid, and viruses were observed using electron 
microscopy gEM-2200FS/CR). 

TCID50 titration and characterization of isolates HY12 

Titration of TCID50 for HY12 isolates was performed using 96 
weU-plates. Briefly, viruses were diluted at a lOxserial dilutions 
and used to infect the quadruplicate wells for each dilution. 48 h 
post inoculation, the cytopathic effects were observed and counted, 
and the TCID5Q was calculated following a standard proce- 
dure [30] . Biological and physiochemical properties of H Y 1 2 were 
characterized according to the standard protocols. HY12 viruses 
were treated either with organic solvents (chloroform and ether), 
heated at 50°C, 60°C, 70°C, and 80°C for 1 h, or incubated at 
pH 3.0 and pH 5.0 before infecting the cells. The stability to heat, 
acid or organic solvents were determined by comparing the 
TCID5Q for treated groups and untreated controls. 

RNA isolation, cDNA synthesis and PCR amplification 

RNA isolation was performed as previously reported[31]. 
Briefly, the infected cells were lysed in TRNzol reagent (Tiangen, 
Beijing), then mixed with 0.2 volume of chloroform and shaken 
vigorously for 30 sec. After centrifugation at 1 2,000 xg for 20 min 
at 4°C, the aqueous phase was mixed with equal volume of 
isopropanol, centrifuged at 1 2,000 xg for 20 min at 4°C. The 
pellet was washed with 70% ethanol and dissolved in DEPC- 



Table 2. Primers used for amplifying the complete genome sequence for HY12. 





Fragment 


primer sequence (5'-3') 


positions 


F1 


S TTAAAACAGCCTGGGGGTTGTA 


1-23 




AS GGTGAAATTTGGTAGCATTGCACT 


1379-1356 


F2 


S CAGTGACACCGACGCAACATC 


1150-1170 




AS CTAAAGTACATAAGCCCGAGAAAGT 


3159-3135 


F3 


S TGACGAGAGCATGATCGAGAC 


2644-2664 




AS CATCTCAGTGAATTTCTTCATCCA 


4011-3988 


F4 


S AATGGGTTCCGATTCCGTTGTG 


3886-3908 




AS AGTCAATTCCAGGGAGGTATCAG 


5632-5610 


F5 


A TCGCTGTCATTCAGAGTGTTTCC 


5262-5284 




AS GACGGCTTTTTCCTTTCTTGATCTT 


6485-6461 


F6 


A AATGTGCAAAGAGACATGCCAGAG 


6455-6178 




AS ACACCCCATCCGGTGGGTGTAT 


7416-7395 



S stands for sense; AS refers to antisense. 
doi:l 0.1 371/journal.pone.0097730.t002 
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Figure 1. Cytopathic effect of HY12 in ceil culture and EM observation. HY12 caused a typical cytopathic effect in MDBK cells after 6-8 h 
inoculation. Cells became rounded with an increased refraction (B). 24—48 h post infection, majority of the infected cells detached off the flask (C). 
The normal MDBK cells were used as negative controls (A). HY12 virus particles observed by electron microscopy to be about 22—28 nm in diameter 
as indicated by arrow (D), the scale bar is 100 nm. 
doi:1 0.1 371 /journal.pone.0097730.g001 



treated HjO. The resulting RNA was kept at — 80°C for further 
analysis. 

The reverse transcriptase reaction was performed using 
Superscript II Reverse Transcriptase (Invitrogen, Carlsbad, 



CA). Briefly, cDNA was synthesized in a volume of 20 H-l 
containing 25 mM Tris-HCI, pH 8.3, 37.5 mM KCI, 1.5 mM 
MgCI2, 5 mM DTT, 0.25 mM each of dATP, dCTP, dGTP and 
dTTP, 40 units of RNase inhibitor, 200 units of M-MLV reverse 



Table 3. Partial physiochemical properties for HY12 strain. 





Treatment 




treated groups (TCID50/0.01 ml) 


untreated groups (TCID50/0.01 ml) 


chloroform 






10-'°' 


ether 




10-'°° 


,0-10.2 


Heat-resistant Testd h) 


50°C 


10-=" 


,0-10.7 




60°C 


10-°° 


,0-10.0 




70°C 


10-°° 


,0-10.4 




80°C 


10-°° 


,0-10.5 


Acid-resistant(pH) 


5.0 




,0-9.6 




3.0 


10-'-' 


,0-10.4 



doi:l 0.1 371 /journal.pone.0097730.t003 
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Figure 2. Unique amino acid mutations in the capsid proteins encoded by I-IY12 isolates. The amino acid sequence of I-IY12 isolates were 
deduced from the nucleotide sequence, and was aligned with 12 known BEV strains in the GenBank. Alignment analysis was performed using each 
HY12-encoded structural protein as a template. Results were shown respectively for VP1 (A), VP2 (B), VP3 (C), and VP4 (D). The identical amino acids 
were marked with symbol and different amino acids to HY12 were presented as individual amino acid symbol. The unique mutation for HY12 was 
highlighted with red color. Deletion of amino acids were marked as 
doi:1 0.1 371/journal.pone.0097730.g002 



transcriptase, 2 |Xg of total RNA, and 2.5 |a,M random primers. 
The cDNA synthesis was carried out at 42 °C for 1 h. 

PGR amplification was done using Taq DNA polymerase 
(Takara Bio Group). The reaction was performed in a total volume 
of 50 111 containing 20 mM Tris-HCI, pH 8.4, 50 mM KGI, 
3 mM MgCIg, 0.25 mM each of dATP, dCTP, dGTP and dTTP, 
5 unit of Taq DNA polymerase, 1 |a.M of each primer, and 2 \xl of 
cDNA synthesized above. The amplification was done after 
optimizing the condition. The primers used for amplifying the 
potential virus sequences were listed in table 1; and primers to 
amplify the complete genome sequence of HY12 isolates were 
listed in table 2. 

Cloning and sequencing 

PGR products were analyzed by electrophoresis using 1% 
agarose gel and cloned to pGEM-T vector (Promega, Madison, 
WI). Recombinants were confirmed by sequencing (Sangon 
Biotech, Shanghai). The resulting sequences were analyzed using 
DNAstar Lasergene software. The nucleotide sequence served as a 
template for searching homologous sequences through GenBank 
(www.ncbi.nlm.nih.gov). 

Alignment and phylogenetic analysis 

Alignment analysis of multiple sequences was performed using 
the Clustal W method [32]. The amino acid sequence of HY12 
was deduced from the nucleotide sequence. Briefly, the nucleotide 
sequences of the 5'-UTR, VPl, VP2, VP3 and VP4 were aligned 
with the corresponding regions of known BEV strains in the 
GenBank, and phylogenetic analysis was performed by neighbor- 
joining methods [33]. 

Results 

Virus isolation and EM observation 

To isolate the potential etiological agents, the faecal samples 
were processed and used to inoculate the MDBK cells, and for 
electron microscopy observation. After incubating with the 
inoculum, MDBK cells showed a typical cytopathic effect as early 
as 6-8 h. Initially, cells became rounded with an increasing 
refraction. 24~48 h post inoculation, majority of the infected cells 
detached off the flask (Fig IB-C). To rule out the possibility of 
toxin effect from the sample, the cultures with inoculum were 
blindly passaged at least 5 generations, and a similar cytopathic 
effect was observed for each passage, indicating the CPE is the 
result of pathogen in the inoculum. To determine the potential 
viruses, the infected cells were frozen and thawed three times and 
processed for electron microscopy observation. As shown in 
Fig ID, virus particles of 22~28 nm in diameter were observed in 
infected ceUs, which was similar to the observation in the faecal 
sample (not shown). The isolated virus was designated as FIY12. 

Virus titration and partial physiochemical properties 

To determine the infectivity of HY12, TCIDjo was determined 
as previously described [30]. Two days post inoculation, wells 
(cells) with cytopathic effect were counted. Results from three 
repeats showed the TGID50 for HY12 isolates is 10 ""^''/O.l ml. 
To further characterize the HY12, partial physico-chemical 



properties, heat-resistant, and acid-resistant experiments were 
performed. As shown in table 3, TCID50 for HY12 isolates has no 
significant change after treating with either chloroform or ether, 
suggesting it is a non-enveloped virus. The TCID50 has no 
significant change for treatment at 50°C for 1 h; however viruses 
completely lost infectivity at 60°C for 1 h, indicating it is sensitive 
to heat treatment at 60°G. As shown in table 3, HY12 isolates is 
stable to acid treatment at pH 3.0. Taken together, the above 
results suggest that HY12 is likely a picomavirus. 

Molecular characterization and the complete genome 
sequence of HY12 

To characterize the FTY 12 isolates, the primers for bovine 
enteroviruses were designed and used to amplify the potential virus 
genome sequence. Simultaneously, the primers for bovine 
parvovirus (BPV) and foot and mouth disease virus (FMDV) were 
used as negative controls. As expected, no fragments for BPV and 
FMDV were amplified after PGR amplification. However, a 
fragment with expected size was obtained with the primers for 
enteroviruses. Further cloning and sequencing this fragment 
turned out that the fragment consists of a sequence with high 
homology to a bovine enterovirus strain SL305, confirming 
FrY12's status as an enterovirus. 

To further characterize the \irus, the complete nucleotide 
sequence of HY12 isolates was dett-rmined using several pair 
primers (table 2). After sequmicing and assembling the PCR- 
ampKfied overlapping fragments, the complete genome sequence 
of HY12 was revealed to consist of 7416 nucleotides, with a typical 
picornavirus genome organization including a 5'-UTR, a large 
single ORE, and a 3'-UTR. The ORF is located between 
nucleotides 818 and 7348, encoding a polyprotein of 2176 amino 
acids with a predicted molecular weight of 243 kD. Comparison of 
the nucleotide sequence revealed that HY12 strain has a similar 
length of nucleotide sequence for 5'-UTR and 3'-UTR to known 
BEV strains. The complete nucleotide sequence of FTY 12 isolates 
was deposited in GenBank (KF748290). 

Unique amino acid mutations in the capsid protein 
encoded by HY12 

To analyze the protein encoded by HY12, the deduced amino 
acid sequence of HY12 was aligned with additional 12 BEV 
strains. Analysis of deduced amino acid sequence for FTY 12 
revealed a few highly conserved regions including nonstructural 
proteins 2A, 2B, and P3 in relation to other bovine enterovirus 
(not shown). Similar to those observed in other enterovirus E, the 
structural proteins encoded by HY12 contains several relatively 
conserved regions (Fig. 2). In addition, there are also several 
variable regions and unique mutations in the capsid proteins, 
especially in VPl and VP4. The mutations include an insertion of 
5 aa at the N-terminal of VP 1 ; an insertion of 2 aa and a deletion 
of 1 aa at the middle region of VTl; a unique insertion of 2 aa (NF) 
at G-terminal region of VPl; and a deletion of 1 aa in the G- 
terminal for VP3 in relation to the known enterovirus E or F. 
Furthermore, aligrraient analysis also revealed rare point muta- 
tions within the VP4 for FTY 12 strain in relation to knovra 
enteroviruses E or F (Fig 2). 
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Figure 3. Pliylogenetic analysis clustered HY12 strain to a new serotype/genotype within enterovirus E. Phylogenetic tree were 
generated by neighbor-joining methods by comparing the sequence regions of 5'-UTR, VP1, VP2, VP3, VP4, 3D, and the complete genome sequences 
for 15 enteroviruses. HY12 strain was placed to the cluster of enteroviruses E after phylogenetic analysis with the all nucleotide sequence regions 
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except 5'-UTR (B-F). The HY12 strain was revealed as a new serotype/genotype (serotype/genotype 3) that only consists of D1 4/3/96 and HY12 
strains in relation to serotype 1 (LC-R4, VG5-27,Vir 404/03) and serotype 2 (SL305, K2577,PS 42, PS 83) enterovirus strains (B-E). When nucleotide 
sequences for the non-structural proteins 3D and the complete genome sequence were employed, the HY12 were clustered to the same clade most 
close to SL305 and K2577 within enteroviruses E (F). However, HY12 strain was clustered to neither clade in enteroviruses E nor enterovirus F using 
the 5'-UTR sequence (A), suggesting an intraserotypic recombination during HY12 evolution. The position of HY12 was highlighted with a triangle. 
doi:1 0.1 371/journal.pone.0097730.g003 
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Figure 4. Recombination revealed In the HY12 strain. Neighbor-joining trees of the structural proteins VP1-VP4, and the non-structural protein 
3D of 15 enteroviruses were compared. When amino acid sequences for VPl, VPS, and VP4 were used to generate phyiogenetic tree, similar patterns 
were observed as those in Fig 3B, 3D, and 3E, indicating a interserotype recombination for the HY1 2 strain. Like the observation in Fig 3A, F and G, the 
HY12 was clustered closely to K2577, SL305, PS 42 and PS 83 strains, an indication of complex interserotypic and intraserotypic recombination in the 
evolution for HY12. The positions of HY12 were highlighted with a triangle. 
doi:1 0.1 371/journal.pone.0097730.g004 
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HY1 2 belongs to a novel serotype/genotype in bovine 
enteroviruses 

To define the relationship of HY12 with other enteroviruses, the 
nucleotide sequences of the 5'-UTR, VPl, VP2, VPS, VP4, 3D, 
the complete genome sequence of HY12 were aligned with the 
corresponding regions of known BEV strains, and phylogenetic 
trees were generated by neighbor-joining methods. As shown in 
Fig 3, when the nucleotide sequences encoding the structural 
proteins (yPl,VP2, VP3,VP4) for HY12 were used to generate the 
phylogenetic tree, the HY12 strain was clustered to entercniruses 

E, but was neither in serotype/genotype 1 (LC-R4, VG5-27,Vir 
404/03) nor serotype/genotype 2 (SL305, K2577,PS 42, PS 83), it 
belongs to a new serotype/genotype consisting of the D14/3/96 
and HY12 (Fig 3B-E). We named it as serotype/genotype 3. It is 
interesting to note when nucleotide sequences for the non- 
structural proteins (3D) and the complete genome sequence were 
employed, the HY12 were clustered to the same clade most close 
to SL305 and K2577 within enteroviruses E (Fig 3A, F, and G). 
However, when the 5'-UTR sequence was aligned with other 
enteroviruses, the sequence identity of FrY12 was 71.7~86.7% 
with other enteroviruses, and phylogenetic analysis showed the 
HY 1 2 is neither in the clade for enteroviruses E nor in enterovirus 

F, it is in a clade with SL305, close to K2577, another separate 
clade (Fig 3A). As shown in Fig 4, similar expected patterns were 
revealed when amino acid sequences for VPl , VP3, and VP4 were 
used to generate phylogenetic tree. Like the obsen'ation in Fig 3A, 
F and G, the HY12 vvas clustered closely to K2577, SL305, PS 42 
and PS 83 (Fig 4B). The above r(-sults indicated a complex 
interserotypic and intraserotypic recombination in the evolution 
for HY12. 

Discussion 

In this study, we reported the isolation of an enterovirus ITY12 
strain from a cattle farm with an unknown disease of high 
morbidity and mortality. Sequencing results demonstrated that 
HY12 is a new isolates within enterovirus E. As far as we know, 
this is the first report of enterovirus E in China. Enteroviruses E 
and F both belong to the genus enterovirus that are etiologically 
associated with bovine enterovirus infections with clinical signs 
varying from respiratory diseases to enteritic, reproductive disease 
and infertility [2,3,4,5,6,7,8,9,10,11,12]. .Although there are 
reports neglecting the pathogenicity of BEVs, increasing reports 
for isolation of BEVs from respiratory and diarrheic cattie suggest 
the relatedness of BEV with those diseases. Our detection and 
isolation of bovine enterovirus FrY12 in a cattle farm with severe 
respiratory and enteric diseases clearly indicate that HY12 might 
play a critical role in this outbreak as etiological agents. 

The classification of bovine enteroviruses undergoes an 
evolution since their initial discovery. Zell reported a molecular- 
based BEV classification system based on virus genetic variability 
and molecular difference, and demonstrated a congruence of this 
molecular-based classification with previous classification methods 
[29]. Based on the generally accepted definition for picomavirus 
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species and serotype, BEVs were classified into enterovirus E 
(BEV-A) and enterovirus F (BEV-B) by comparing sequences from 
5'-UTR, 3'UTR and capsid protein region, where a different 
serotypes were further divided for either enterovirus E or 
enterovirus F [27,28,29]. It is generally accepted that picornavirus 
serotypes are molecularly defined by tlu; diversity of the capsid 
proteins, while the enterovirus species were determined by the less 
diverse non-structural protein regions. The percentage of sequence 
identity was used for species/ serotypes definition, where a range 
from 50% to 55% for heterologous species, 70% to 85 % for 
heterologous serotypes/homologous species, and greater than 
90 % for homologous serotypes [29] . According to this criteria, we 
defined FrY12 enterovirus isolates as enterovirus E since the amino 
acid sequence identity of VPl encoded by HY12 is only ranged 
from 53'X) to 55% with previously identified BEV-B viruses such as 
BEV strains PS87, 3A, 261, and two newly strains isolated in 
China (BHM26, and BJOOl). The finding ti^at VPl of HY12 is 
about 80-83% with those previously defined BEV-A suggests 
HY12 is a heterogeneous serotype within enteroviruses E. 
Moreover, phylogenetic analysis clearly showed that HY12 isolates 
belongs to neither serotype/genotype 1 (SL305, K2577 PS42, 
PS83) nor serotype/genotype 2 (LC-R4, Vir 404-03 PA12-24791); 
it belongs to a new serotype/genotype. The results that HY12 is 
phylogenetically clustered with D 14/3/96 strain [29], a strain 
reported to be difficult in serotype/genotype typing further 
support the classification of FrY12 to new serotype/genotype 
within enterovirus E. 

Interserotypic and intraserotypic recombination have been 
defined previously in poliovirus, echovirus and enteroviruses 
[1,34,35,36,37]. The incongruence between phylogenies of 
different genome regions is considered as an indication of 
recombination events in enteroviruses [29,38]. Phylogenetic 
analyses of the 5'-UTR, VPl, VP2, VP3, VP4, and 3D of 
HY12 demonstrated the incongruence between the phylogenies, 
indicating a complex interserotypic and intraserotypic recombi- 
nation in the evolution of HY12 isolates. 

In conclusion, this study reports the isolation and characteriza- 
tion of a novel enterovirus HYl 2 strain from a severe outbreak 
characterized with respiratory and enteric disease in a cattie farm 
with high mortality and morbidity in China, and provides the 
molecular evidence for defining HYl 2 strain as a new serotype/ 
genotype 3 within enterovirus E. 
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